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A general expression of the current induced spin torque in a magnetic layered structure in the 
presence of external dc or ac voltages is derived in the framework of the scattering matrix approach. 
A detailed analysis is performed for a magnetic-nonmagnetic-magnetic trilayer connected to external 
leads in the presence of dc voltage bias in the ballistic regime. Alternatively, the possibility of 
producing spin torque by means of the adiabatic ac modulation of external gate voltages (quantum 
pumping) is proposed and discussed. 



PACS numbers: 73.23.-b,72.25.Pn,75.60.Jk,72.15.Qm 



INTRODUCTION 



Multilayers of alternating magnetic (generally ferro- 
magnetic) and nonmagnetic metal layers have recently 
attracted a lot of attention because of giant magnetore- 
sistance (GMR) effects. In fact their electrical resistance 
depends strongly on whether the moments of adjacent 
magnetic layers are parallel or antiparallel and this ef- 
fect has allowed the development of new kinds of mag- 
netic memory devices ^ The origin of the GMR effect is 
the stronger scattering of conduction electrons by a mag- 
netic layer when their spins lie antiparallel to the layers 
magnetic moment compared to the case when their spins 
are parallel to the moment. Thus the orientations of 
magnetic moments can affect the flow of electrons, but 
as a reciprocal effect a polarized electron current scat- 
tering from a magnetic layer can affect the moment of 
the layer itself. In fact, as proposed by Berger— and 
Slonczewski, - an electric current passing perpendicularly 
through a magnetic multilayer may exert a torque on the 
moments of the magnetic layers. This effect is known as 
"spin transfer, and may alter the magnetization state of 
the layer. Of course this is a different mechanism from 
the effects of current induced magnetic fields and offers 
the possibility of realizing new kinds of magnetic devices. 
Importantly, it can serve as a mean to detect the spin 
dynamics at nanometer scale and measure spin currents. 
However, in order to utilize these effects in real devices, 
it is necessary to achieve a quantitative understanding of 
spin-current-induced torques and also provide different 
ways for their generation. 

In the present literature a semiclassical WKB approx- 
imation with spin dependent potentials in ID has been 
proposed in Ref.[3]. A generalization of this study has 
appeared in Ref.[4] and an extension of this calculation 
to take into account the band structure effects on the 
degree to which an electron is transmitted through a 
magnetic/nonmagnetic interface has been proposed by 
Brataas et al£ by means of the kinetic equations for spin 
currents and by Waintal ct al. 6 by using a Landauer- 
Biittiker type of approach. A first principle theory based 
on the spin mixing conductance appeared in Q ■ 



Here we present a general formalism based on the scat- 
tering matrix approach to calculate the torques in a mag- 
netic layered structure in the ballistic regime. Then we 
use this formalism to make an explicit calculation for 
a magnetic(M)-nonmagnetic(NM)-magnetic(M) trilayer 
connected to two metallic external leads. In this case 
the ballistic regime is characterized by a spin-diffusion 
length l s and a mean-free-path l m larger than the NM 
spacer. We perform first the calculation in the presence 
of a dc voltage bias applied to the leads, or equivalently 
for currents perpendicular to the microstructure, then we 
proposed an alternative way of achieving a spin-torque by 
the adiabatic ac modulation of two system parameters (in 
our case the barrier heights at the M-NM interface), i.e. 
by a quantum pump mechanism 8 . Although we limit our 
calculation to the single channel case, differently from 
Refjy all spin-scattering processes at the M-NM inter- 
face are properly taken into account, including spin- flip 
mechanisms. 

The organization of the paper is the following: In SeclHl 
we introduce the model Hamiltonian and the general ex- 
pression for the spin-torque. In Sec lIIII we develop the 
scattering matrix approach generalized for the calcula- 
tion of the spin-torque in a magnetic layered structure. 
In Sec|IV]we discuss the adiabatic quantum pumping of 
spin-torque and finally in Sec|V] we present the results 
for a specific system, a trylayer of M-NM-M metals con- 
nected to two nonmagnetic leads. A quantitative analysis 
of the spin-torque components is presented both in the 
case of a dc bias applied to the external leads and when 
a quantum pump mechanism is activated. The two ways 
of generating a spin-torque presents some differences that 
are comparatively discussed. These differences could be 
important for the realization of new magnetic devices. 



II. 



THE MODEL AND FORMALISM 



Let us consider a ID system whose Hamiltonian is: 

H = ~im®x + ^f 1 !^! where V(x,a) is some spin- 
dependent potential operator. The field operator describ- 
ing the electron state 9, obeys the Schrodinger equation 
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= ihd t ^ and = -ihd t ^ and the time evolu- 

tion of the electron charge density p — is deter- 

mined by the Heisenberg equation: 



d t (^eV) = j[&(H- H)V] = -d x J, (1) 



where J = ^^ev^ is the current density, while v — 

■^r(9 x — d x ) is the velocity operator along the x- 
direction. From ([1]) the well known continuity equation 
dtp + d x J = is obtained. 

In an analogous manner, the time evolution of the spin 
density = ^ja^, p = {x,y,z} (a^ are the Pauli 
matrices) , can be derived: 

d t Sa = ~d x J s u + Tn. 



(2) 



Here J* = ^^a^v^ is the /j-th component of the spin 
current density and is the the density of spin torque. 
When the spin-dependent potential describes the Zecman 
interaction between the spin and a magnetic field along 
the direction n{x), V(x,d) = r y(x)h(x) ■ er, the density 
of spin torque is explicitly given by 7), = ^{x)^ (n{x) x 

= ^T-(H X ) x SV When >y(x) = 0, the spin 
torque term vanishes and a continuity equation is obeyed 
by the spin currents, d t S + d x J s = 0. In the case of a 
system in which a magnetic central region x £ [—a, a] is 
connected to two external non-magnetic leads, under the 
stationary condition for the spin-current (i.e. dtS = 0) 
the spin torque per unit of area can be calculated by 
spatial integration of the equation d x J* = r M along the 
x-direction: 



a = {left, right} represents the lead index, while the 
S- matrix must be unitary, i.e. S^S — 1. The quantum 
field representing the electron state in the lead a can be 
expressed in terms of the scattering operators as follows: 



E 



[e lkx a a a {E)+e 



dEp a (E) exp 
lkx b a a {E)l 



E 



(5) 



where p a (E) = [y/2Trhv a (E)] 1 is the density of states 
of lead a, v a (E) is the velocity of the electrons with wave 
vector k — \ / 2mE/h and the velocity is chosen to have 
positive orientation for the incoming states. The spin- 
current density J s = ^ ^av^ can be calculated with 
the approach of Ref.~ and the result is: 



js = J_ 

4tt 



dEdE' exp 



(6) 



[a a \E)a^a a {E') - b a \E)a^b a {E% 



where the following spinorial representation has been in- 
troduced: 



while 



/at at\ 



(7) 



Using the relation 



{a^{E)4,{E')) = Saps,?,?' fp(E), f (E) being the Fermi 
function of the lead (3 and Eq.Q, the p-th component 
of the spin density current ( J* a ) in the lead a can be 
written as^: 



7 s — 7 s 



+ a 



dxf = T, 



(3) «J = 



where J S R / L — J s (±oo) is the value of the current at in- 
finity. The above result tells that the spin torque per unit 
of area is due the difference of the spin currents on the 
right and on the left of the magnetic scattering region. In 
the following Section we analyze the spin torque transfer 
within a scattering matrix approach. 



III. DENSITY OF SPIN CURRENT WITHIN 
THE SCATTERING APPROACH 

In order to calculate the spin torque through a mag- 
netic region we employ the scattering matrix formalism 
of Ref.[10( and consider a single channel approximation 
to simplify the derivation. 

The S'-matrix connects the outgoing states to the in- 
coming states through the relation: 



E*; 



aa' u a' 



(4) 



where b% , a" are the scattering operators for the out- 
going and incoming states, a = ± is the spin index, 



4tt 



E 





dETr{S a ^ (E)a„ S a ? {E)}f p (E) ,(8) 



while using ( \5§ the spin torque per unit of area T p = 
/ dxT^ is: 

T, = ^ E [ dETr{S^(E)a,S^(E)}f (E).(9) 

When a dc voltage bias V is applied to the external leads, 
changing their chemical potentials in pi >r = Ep ± eV/2, 
the variation of the spin-torque with respect to the ex- 
ternal perturbation w — eV/2 can be written as S W T^ — 
d w T^5w in the linear response regime. It is determined 
by the torkance d w T^. Using Eq.@ and taking the zero 
temperature limit, we explicitly have--: 



8 T 



1 

2^ 



Tr{a^S 21 S 21 ^-S 12 S 12 ^} 



E — Ef 



(10) 



The torkance along the direction of the unit vector n — 
(n x , n y ,n z ) is obtained as n ■ d w T. Eq.fTU]) describes the 
magnetic answer of the system to an external dc voltage. 
Let us note that compared to previous works^i our calcu- 
lation of the scattering matrix includes also off-diagonal 
spin- flip reflection and transmission amplitudes. 
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IV. ALL ELECTRICAL GENERATION OF SPIN 
TORQUE BY QUANTUM PUMPING 

In this section we discuss the generation of spin torque 
by means of ac external gate voltages. In particular we 
generalize the notion of quantum pumping of charges due 
to Thouless^ to the spin torque. In a quantum pump a dc 
particle current is generated by the ac adiabatic modula- 
tion of at least two out-of-phase independent parameters 
of the system (e.g. local magnetic fields or gate voltages) 
in absence of bias. In our calculation we will show that 
since pumping procedure in a magnetic layered structure 
can generate spin currents other than charge currents, a 
spin-torque is generated by the gradient of spin current 
(see Eqs.([2|)-([3])). In particular, we focus on the system 
of Fig.JT]) in which a microstructure made of a central 
non-magnetic-region (NM) is connected to two external 
non-magnetic leads (LL and RL) through magnetic lay- 
ers (Ml and M2) (whose width is taken less than the De 
Broglie wavelength). Applying the idea of pumping, we 
modulate in time the barriers heights at the interface be- 
tween the M-NM regions by the top gates Gl and G2. In 
the presence of this ac modulation the scattering matrix 
S depends explicitly on time and the relation between 
incoming and outgoing states is: 



b a (t)=J2 J dt'S al3 (t,t')aP(t'), 



(11) 



where the spin-indices are absorbed in the spinorial no- 
tation (see Eq.([7])). When the gates are varied adia- 
batically in time, an instantaneous approximation can 
be made, i.e. S af3 {t,t') = 5{t - t')S af3 {t), and thus 
b a {t) = S aP (t)aP(t). In particular, since the time 
dependence of the scattering matrix S(t) is induced by 
two external parameters of the form 



X 1 (t) 
X 2 (t) 



X% + X? sin(urf) 
X% + X% sin(ut + (p), 



(12) 



then S(t) = S(Xi(t),X2(t)). In particular, when the 
amplitude of the ac parameters is small, i.e. Xf 2 -C X® 2 , 
the scattering matrix can be expanded as follows: 

S af3 (t) ~Sf+ sfe 1 ^, (13) 
i)=±i 



where u> = 2itv is the angular frequency of the adiabatic 

are given by 



modulation and the matrices 



IT) 

~~2 



X?(d Xl S)o + XZe l ^(dx 2 S) 
The Fourier transform of ([TBI is then: 



(14) 



S aP {E) = 2Ti[Sf8{E) 



■q=±l 



sf5(E + r,cj)}. (15) 



Consequently, the relation between the outgoing and in- 
coming states in the Fourier space takes the following 
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FIG. 1: Schematic representation of the NM/M/NM/M/NM 
system described in the main text. 



form: 

b<*(E) = £ [sfa^E) + J2 *fa (E + w)] • (16) 
/9 n 



Using Eq. (|16|) in (J6J) the /i-component of the spin torque 
per unit of area is given by: 

T ^ ^ E / dETr{sf^„sf}f (E + w ). (17) 

r/a/3 

Since no bias is present between the leads fp(E) — f(E), 
and in the zero temperature limit, the /u-th component of 
spin-torque per unit area to leading order in the adiabatic 
frequency u> is: 



(18) 



(all terms independent from the external perturbation 
have been dropped). Using (I14p this can be rewritten 
in terms of the parametric derivatives of the scattering 
matrix as : 



8^ 



J2Tr{Af + Af% (19) 

a(3 



where we introduced the quantity A" 13 = 
i(dx 2 S a ^) a fi {d Xl S a(3 ) . The torque generated in 
the direction of the unit vector n is given by Tn = n ■ T. 
Eq. (fl9|) represents the torque pumped in a magnetic 
system and in principle could be very different from the 
one induced by a dc voltage. In the following we perform 
the explicit calculation for the system of Fig. 1. 

V. SPIN TORQUE IN M— NM— M SYSTEMS 

Here we analyze the spin torque generated in the sys- 
tem of Fig. 1 in the following two situations: (i) in the 
presence of a dc voltage applied to the external leads; 
(ii) in the absence of bias and using a quantum pumping 
procedure. The central NM region is connected to two 
external leads via the thin magnetic layers Ml and M2 
and without loss of generality we choose the magnetiza- 
tion of M2 parallel to the z-axis, while the magnetization 
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of the region Ml has arbitrary direction. The distance 
between the two magnetic barriers is taken 2a and their 
height is controlled by the the top gates Gl and G2. A 
minimal model for the system above is given by the fol- 
lowing Hamiltonian: 



H 



2m 



d 2 x + V m (x,a) + U(x), 



(20) 



where U(x) = giS(x + a) + g2^(x — a) is the barriers 
potential at the interface between M-NM region, while 
V m (x, (?) = 7;<5(x + a)n ■ a + j r S(x — a)a z represents the 
magnetic interaction between the conduction electrons 
spin and the magnetic moment of the layers Ml and M2. 
The scattering problem for the Hamiltonian above can 
be easily solved by imposing the appropriate boundary 
conditions on the wavefunctions and their derivatives at 
the interfaces (see, for instance, Ref . jl2| ) and from the 
knowledge of the scattering matrix the spin torque T M can 
be calculated using the equations of the previous sections. 
The following consideration is in order here: The pres- 
ence of two magnetic layers has the effect of allowing for 
multiple scattering of the electrons between them, which 
gives rise to an explicit asymmetry of the spin current 
and produces a finite spin-torque effect. 

In the following we employ the adimensional quanti- 

2mql Iraqi i \ 2m->v T" 1 2m7j 

ties: z i = Z2 = 7&t> Tr = ft^fc?' Tl = 7F57> 

where kp is the Fermi wavelength, m — am e is the effec- 
tive electron mass, while a represents the ratio between 
the effective and the bare electron mass m e . The dis- 
tances are made adimensional by multiplying by kp, i.e. 
a — > kpa. 



A. Spin torque by means of dc voltage bias 

Here we present the results of the spin torque per 
unit of area generated by a small dc voltage bias V ap- 
plied to the external leads. In the following we use po- 
lar coordinates for the magnetic moment Mi — Min, 
h = (8111(6*) cos(<^>), sin(6 l ) sm(0), cos(0)), and assume that 
the momenta Mi and M% lie in the plane x — z. In this 
case we have n — (sin(0), 0, cos(0)), where the 8 repre- 
sents the angle between Mi and M^. 
In Fig- (J2J) the three components of the torque (in unit 
of eV/2) are shown as a function of angle 8 between the 
magnetic momenta of the regions Ml and M2 and by fix- 
ing the remaining parameters as follows: a — 2, Z\ = 0.5, 

22 = 1,1^ = 0.1,1^ = 1. 

Similarly to other works^, we find that all the torque 
components vanish when Ml and M2 are parallel or anti- 
parallel. Furthermore, the z-component T z of the torque 
is very small compared to T x and T y , while the compo- 
nent perpendicular to the plane of the magnetic moment 
T y is the strongest. This component of the torque fol- 
lows a sin 8 behavior for the specific values of magnetic 
interactions considered in the figure (Tr ~^>Yl). In fact 
when the non-spin-polarized electrons are incident on the 
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FIG. 2: Torque T x (dashed line), T y (dashed-dotted line), T z 
(full line) in unit of eV/2 computed as a function of the angle 
between Ml and M2. The parameters have been fixed as 
follows: a = 2, z x = 0.5, z 2 = 1, T L = 0.1, T R = 1. When the 
magnetization directions of the layers Ml and M2 are parallel 
or anti-parallel the spin torque vanishes. 



magnetic layer M2, spin- filtering removes the component 
of the spin angular momentum perpendicular to the layer 
moment from the current acting as a polarizer along the 
z direction. Thus the polarized electrons scattered off 
M2 and incident on Ml generate an effective torque on 
the layer momentum Ml proportional to the sin 8. The 
negative sign of the spin torque component T y is due to 
the fact that electrons with spins parallel to the moment 
of the magnetic region M2 have a larger transmission 
probability compared to those antiparallel. 

In Fig. ^ the components of the torque are shown as a 
function of the semi-distance a between the gates Gl and 
G2 and by fixing the remaining parameters as follows: 
8 = 1.5, 2i = 0.5, z 2 = 1,T L = 0.1, T R = 1. For a meso- 
scopic system characterized by a De Broglie wavelength 
A ~ 20 — 30nm the adimensional distance a = 1 would 
correspond to 3.18 — 4.78nm. As shown in the figure, the 
torque presents a characteristic oscillatory behavior with 
the semi-length of the central region a. These oscillations 
can be regarded as a quantum-size effect. They reflect 
the perfect ballistic regime of electron transport through 
the spin valve. The physical mechanism behind the os- 
cillations is the interference effect of the electrons propa- 
gating across the non-magnetic-magnetic interface from 
the right lead to the left lead and electrons propagating 
backwards. The particular value of the oscillation follows 
from the values of the spin-dependent Fermi wavevector. 

In Fig. ([J) we plot T x<y>z as a function of the strength 
Z2 of the right barrier (controllable via the gate G2) by 
fixing the other parameters as: 8 — 1.5, z\ — 0.5, a = 2, 
Tl = 0.1, = 1. The torque of course depends cru- 
cially on the transparency of the barrier and thus the 
components T x ^ y decreases with increasing barrier height 
becoming comparable to T z . An interesting feature ap- 
pears by increasing the Zeeman interaction as shown 
in Fig.©: for values of above 0.2 the ^-component 
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FIG. 3: Torque T x (dashed line), T y (dashed-dotted line), T z 
(full line) in unit of eV/2 computed as a function of the semi- 
distance a between the gates Gl and G2. The parameters 
have been fixed as follows: 8 — 1.5, Z\ = 0.5, zi — 1, Tl = 
0.1, Tij = 1. Notice that the torque T y , i.e. the torque in 
the direction perpendicular to the plane x — z presents large 
oscillations compared with the in-plane components. 
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FIG. 4: Torque T x (dashed line), T y (dashed-dotted line), 
T z (full line) in unit of eV/2 computed as a function of the 
non-dimensional barrier strength Z2 controlled by G2. The 
parameters have been fixed as follows: 9 = 1.5, zi = 0.5, 
a = 2, F L = 0.1, T R = 1. 



of the torque T z starts to assume relevant values. In fact 
the magnetic interaction becomes more effective in align- 
ing the electron spins along the magnetic moment of the 
layer Ml thus reducing the contribution to the gradient 
of spin current. 



B. Spin torque generation by means of adiabatic 
quantum pumping 

An alternative way to generate a spin torque is the adi- 
abatic quantum pumping technique. As explained above 
we keep the two external leads at the same chemical po- 
tential and in order to generate a spin current we modu- 
late out-of-phase in time the two barrier heights z\ and 



FIG. 5: Torque T x (dashed line), T y (dashed-dotted line), 
T z (full line) in unit of eV/2 computed as a function of Tl 
controlled by the magnetic momentum in the region Ml. The 
remaining parameters have been fixed as follows: 6 = 1.5, 
zi = 0.5, z 2 = 0.8, a = 2, T R = 1. 



Zi by the gate voltages Gl and G2: 

1 + sin(wi) 



*i(t) 
z 2 (t) 



(21) 



z 2 + z 2 sin(w£ + ip). 



Thus by using (fl"9|) the fi — th component of the torque 
in the adiabatic regime is given by: 



T, = 



'1 ^2 



sin(</?) 



8tt 



a/3 



Tr{Af + Afi}. (22) 



In the weak pumping regime (i.e. zf 3> zf) considered 
here, presents a sin(</j) behavior with respect to the 
pumping phase ip and thus in the following analysis we 
set ip = 7r/2 and measure the spin torque in unit oiKv/2. 
Concerning the frequencies of the pump, we can safely 
consider as adiabatic the frequencies for which hu> is 
smaller than the first gap in the electron energy spec- 
trum such that the system lies in its ground stated. For 
simplicity by considering the energy spectrum of a par- 
ticle confined in a box, we easily find that the threshold 
frequency is v ~ 2m* l' 2 > wnere L is the length of the 
system, while m* = am e is the effective mass of the elec- 
trons. Thus for a system with L = 1.5/^m and effective 
mass ratio a = 0.023, v w 43.4 GHz which is consistent 
with the values of the frequencies of quantum pumping 
in quantum dots. For a smaller system characterized by 
L = 0.5/im and a = 0.023 we obtain i> w 396.5 GHz 
which is a more convenient value to achieve in experi- 
ments. 

In Fig.© we plot the components as a function of 
the angle 9 between the magnetization Ml and M2 and 
by fixing the other parameters as follows: zf — z^ = 0.2, 
r L = 0.9, z\ = 0.5, z£ = 0.3, a = 1, T R = 1. Com- 
pared to the case of dc bias, the 8 dependence of the 
torque is not a simple sin 9 form because the strength of 
the magnetic interactions Tl and Tr is comparable and 
thus both cos 8 and sin# terms contribute. Moreover as 
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FIG. 6: Torque T x (dashed line), T y (dashed-dotted line), T z 
(full line) in unit of hi// 2 computed as a function of 9 between 
Ml and M2. The remaining parameters have been fixed as 
follows: zf = z% = 0.2, T L = 0.9, z? = 0.5, z\ = 0.3, a = 1, 
F R = 1. 
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FIG. 8: Torque T x (dashed line), T y (dashed-dotted line), T z 
(full line) in unit of hv/2 computed as a function of the static 
part of the non-dimensional barrier strength z\ controlled by 
G2. The remaining parameters have been fixed as follows: 
zt = zt = 0.2, T L = 0.9, z\ = 0.5, a = 1, T R = 1, = 1. 
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FIG. 7: Torque T x (dashed line), T y (dashed-dotted line), 
T z (full line) in unit of hv/2 computed as a function of the 
semi-distance a between the gates Gl and G2. The remaining 
parameters have been fixed as follows: zf = z% = 0.2, r,L = 
0.9, 2? = 0.5, zi = 0.3, F R = 1, = 1. 



in the previous analysis, while the in-plane component 
of the spin torque T x and T z assume very small values 
in the parameter region considered, the T y component 
presents relevant values. Indeed, by setting k, 1 and 
for pumping frequencies v of the order of 10-30 GHz we 
have a torque component T y « 0.05 — 0.15/ieV per unit 
of area, which is of the same order of magnitude ob- 
tained in a similar system under dc bias (see Ref.[l3[). 
In Fig. ([7]) we show the spin torque components T XiViZ as 
a function of the semi-distance a between the gates Gl 
and G2, for the remaining parameters: zf = z% = 0.2, 
r L = 0.9, z\ = 0.5, z\ = 0.3, T R = 1, 9 = 1. Com- 
pared to the case in which the spin-torque is generated 
by an external dc bias, more harmonics are present in 
the torque dependence on a. This behavior can be qual- 
itatively explained by the fact that the torque generated 
by the pumping procedure is related to the parametric 
derivatives of the scattering matrix. In Fig. ([5]) we show 



the spin torque components T x y z as a function of the 
static part of the adimensional barrier strength z\ con- 
trolled by G2, while the remaining parameters have been 



fixed as follows: 



= 0.2, T L - 0.9, z\ 



0.5, 



a = 1, Tr = 1, = 1. As shown, for increasing values of 
the static barrier strength z% the torque pumped in the 
system becomes vanishes. Indeed, the reduction of the 
spin current from the external leads reduces the gradi- 
ent of the current and thus the spin torque is strongly 
suppressed. In the case of pumping, the suppression of 
the torque by increasing the barriers height may be in 
general stronger compared to the dc case. In Fig.© we 
show the spin torque components T x ^ y ^ z as a function 
of the Zeeman interaction Tl for the other parameters: 



0.2, z° 2 



1, z\ = 0.5, a = 1, T 7 



1. 



A behavior of the torque similar to the one presented in 
Fig. ([5]) is found. The magnetic interaction of the electron 
spins with the magnetic moment of Ml causes a suppres- 
sion of the torque above Tl ~ 1-7 thus acting like an 
energy barrier for the electrons. 



VI. CONCLUSIONS 

Within the scattering matrix approach we studied the 
spin torque generated by dc or ac external perturbations 
acting on a multilayered system consisting of a sequence 
of magnetic/nonmagnetic regions. In particular, we stud- 
ied the spin-torque of a magnetic-nonmagnetic-magnetic 
trylayer connected to metallic nonmagnetic leads. We 
have focused on the effect of spin-filtering as the mech- 
anism for current induced torque, i.e. the difference in 
the transmission and reflection probabilities for electrons 
with spins parallel and antiparallel to the moments of 
the magnetic layers. As a source of spin-dependent scat- 
tering we included also the spin-flip mechanism (i.e. the 
off diagonal elements of the scattering matrix) . The spin 
torque generated by the application of a dc bias to the 
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FIG. 9: Torque T x (dashed line), T y (dashed-dotted line), 
T z (full line) in unit of hv/2 computed as a function of Tl 
controlled by the magnetic momentum in the region Ml. The 
remaining parameters have been fixed as follows: zf — — 
0.2, z% = 1, zl = 0.5, a = 1, Vr = 1, = 1. 



external leads was analyzed as a function of the relative 
orientation 9 between the magnetic moments of the mag- 
netic regions, the length 2a of the nonmagnetic spacer, 
the barriers transparencies and the Zeeman couplings. 
From the analysis we observed that the spin-torque van- 
ishes when the magnetic moments of Ml and M2 are 
parallel or antiparallel, while an oscillating behavior of 
the spin-torque was observed as a function of a, indicat- 
ing a quantum size effect. When the transparencies of 
the tunnel barriers at the interfaces between the mag- 



netic/nonmagnetic layers are lowered, a strong suppres- 
sion of the spin torque is observed due to the reduction 
of the spin fluxes through the barriers. Similar effects 
were found as a function of the strength of the magnetic 
interaction T^. 

As an alternative to considering external dc bias, we 
proposed a current induced spin torque based on quan- 
tum pumping, generalizing the original idea of Thouless 
for the charges^. We formulated a scattering matrix for- 
malism of the spin torque pumping. By modulating in 
time the strength of the two out-of-phase tunnel barriers 
at the interface between the magnetic/nonmagnetic lay- 
ers, spin polarized currents were produced in the leads 
and a gradient of spin current responsible for the spin 
torque was generated. By studying the spin torque com- 
ponents in this case we observed different signatures of 
the magnetic reversal mechanism due to the pumping 
procedure. For instance, within the weak pumping limit 
additional harmonics were observed in the oscillations of 
vs the length of the nonmagnetic region 2a and a 
stronger suppression of the spin torque at increasing the 
barrier heights. Thus the pumping mechanism could offer 
an additional tool to study experimentally the magnetic 
response of metallic heterostructures. 

Experimentally, current induced spin torque has been 
realized in Co/Cu/Co sandwich structure o 15 ! 16 and in 
Py/Cu/Py heterostructuresii. Our proposal can also be 
realized by using two ferromagnetic EuS barriers coupled 
to Al leads and separated by an AI2O3 spacer of few tens 
of nanometers by modifying, for instance, the experimen- 
tal setup studied in Ref . [18j . 
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